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The PAMELA, ATIC and Fermi collaborations have recently reported an excess in the cosmic ray 
positron and electron fluxes. These lepton anomalies might be related to cold dark matter (CDM) 
particles annihilating within a nearby dark matter clump. We outline regions of the parameter space 
for both the dark matter subhalo and particle model, where data from the different experiments are 
reproduced. We then confront this interpretation of the data with the results of the cosmological 
N-body simulation Via Lactea II. Having a sizable clump (Vmax ~ 9 kms~^) at a distance of only 
1.2 kpc could explain the PAMELA excess, but such a configuration has a probability of only 0.37 
percent. Reproducing also the ATIC bump would require a very large, nearby subhalo, which is 
extremely unlikely (p ~ 3 x 10"''). It is even less probable for the smaller Fermi bump to be caused 
by the presence of such an object. In either case, we predict Fermi will detect the gamma-ray 
emission from the subhalo. We conclude that under canonical assumptions, the cosmic ray lepton 
anomalies are unlikely to originate from a nearby CDM subhalo. 

PACS numbers: 95.35.-|-d,98.35.Gi,11.30.Pb,95.30.Cq 

Recent measurements [H 0, 0] of cosmic ray positrons and electrons (hereafter dubbed leptons) have stirred a lot of 
interest, since they may be the first indirect hint of the presence of particle Dark Matter (DM) in the halo of the Milky 
Way (MW) 0,i,i01. In the standard cosmic ray picture, positrons are secondary species produced by the spallation 
of cosmic ray protons and helium nuclei on the interstellar medium Q . But secondary positrons quite clearly fail to 
reproduce the PAMELA [ij, ATIC or Fermi Q measurements. Note that H.E.S.S. also measured cosmic ray leptons 
spectrum]^, but with higher threshold, making these measurement less relevant for dark matter interpretations. 
While astrophysical primary positron sources exist: (the most obvious class of candidates being pulsars) and can 
account for the recent measured anomalies [O, 12, Tj], a more exotic origin is possible. In particular, the observed 
excess could be sourced by the annihilations of the massive and weakly interacting particles (WIMPs) proposed to 
explain the astronomical DM. Estimates based on standard assumptions on the annihilation cross-section and on the 
DM halo fail, however, to reproduce the measured cosmic ray lepton anomalies. Enforcing a standard thermal DM 
production in the early universe sets the pair-annihilation rate a couple of orders of magnitude below what is needed 
to explain the data. This issue is usually tackled by assuming boost factors which can arise from different arguments. 

A first possibility lies in the existence of a particle physics boost factor which increases the annihilation rate as a 
result of a mechanism evading the relic density constraint (e.g. velocity-dependent annihilation cross-sections, non- 
thermal primordial production) . Alternatively, an effective astrophysical boost factor could stem from the dumpiness 
of the Galactic DM halo. Indeed, the annihilation rate being proportional to the squared DM density p^, the 
"clumpier" the halo, the higher the signal in a (p^)/ (p)^ proportion. Although the boost factor is actually sensitive 
to the energy [3, 15, [l^, most analyses of the cosmic ray lepton anomalies have so far treated dumpiness either 



assuming an energy-independent enhancement of the signal or modifying its spectral shape regardless of the overall 
normalization. Moreover, the various contributions to the positron flux at the Earth have not been derived withxin 
the same cosmic ray propagation model and are in general inconsistent with each other. 

In this letter, we re-assess 17, 18, [l^ the possibility that a single nearby DM clump contributes substantially 



to the lepton anomalies. Particular attention is paid to cosmic ray propagation. We adjust the clump distance 
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FIG. 1: Best fits to the PAMELA data in the case of a positronic fine (see the e /e 
the ATfC data (right panei). 
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D and luminosity L in order to reproduce the PAMELA, ATIC and Fermi data. The probabilities of these clump 
configurations are calculated based on the cosmological N-body simulation Via Lactea II [2^, which allows us to 
quantify for the first time how unlikely a sufficiently bright nearby CDM subhalo is. We eventually comment about 
the interplay between the DM particle properties and the clump parameters. As an illustration, we point out extreme 
configurations where these mutual effects lead to subtle, but relevant modulations of the DM signal. 

The lepton density "0 = dn/dE is related to its source q through the stationnary cosmic ray diffusion equation 
described in Q . Lepton propagation throughout the diffusive halo is dominated by space diffusion and energy losses 
via synchrotron emission and inverse Compton scattering on the CMB and stellar light. The diffusion equation 
is solved in the framework of the two-zone diffusion model detailed in 2l| and yields the flux (p = ajj/^T^- The 
propagation parameters best-fit the B/C ratio and correspond to model med of [23 • A value of 7 x 10^^ sec is taken 
for the energy loss timescale. 

The positron flux at the Earth 0g+ — (ff^^ + 0^+ -I- (jf^^ results here from three contributions. The astrophysical 
background (jf^+ is provided by the secondary species produced by primary cosmic rays impinging on the interstellar 
material and is computed as in Q . The smooth DM halo contributes the source term 



9dm (x, E) 



Ps(x) 



.f{E) , 



(1) 



where f{E) is the energy spectrum of the positrons created in the annihilation process. The Galactic DM halo density 
Ps is borrowed from the results of the Via Lactea-II simulation, with a spherical profile featuring an inner (outer) 
logarithmic slope of -1.24 (-3) and a 28.1 kpc scale parameter. The local density pq is equal to 0.37 GeV cm""^. The 
galactocentric distance of the Earth is tq = 8.5 kpc. Finally, the contribution of a nearby clump located at can 
be expressed as 



1 



g^M(x,i?) = 77 M <5^(x-x,) f{E) 



(2) 



where L = /^.j^^-jp Pc(x) d'^x is defined as the subhalo luminosity. Furthermore, as clumps are treated here as point- like 
objects [13], the source terms q^^^ and q'^^ add up directly to yield the total DM lepton signal. 

The annihilation cross-section of the DM particles under scrutiny is set equal to {av) = 3 x 10^^® cm'^ s^^. This 
canonical value matches a thermal production of DM in the early universe. We also consider 100 GeV and I TeV 
WIMPs as benchmark cases. Finally, inspired by [1, we concentrate here on leptophilic species and considered 
either a pure e"'"e~ annihilation final state (positronic line) or an equal production of charged leptons e^ -I- -\- r^. 
Notice that a case for both scenarios can be made from the model-building perspective. We also tried to use models 
featuring pure bb or W+W~ annihilation final state, but we disregard them since the required clump configurations 
are extremely unlikely {p < fO~^). 

The particle physics framework being set, we perform fits to the PAMELA, ATIC and Fermi data which include a 
smooth DM component plus a contribution from a DM subhalo whose luminosity L and distance D are free parameters. 
For PAMELA, we compute the positron fraction 0e+ /0e+ -I- (pe- where we use for (p^,- the observed cosmic ray electron 
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flux measured by AMS and HEAT Q and parameterized by Casadei & Bindi As regards ATIC and Fermi, 
we derive the total lepton flux c/)e+ + (f)^,- , assuming that the electron background (j)^'^'^^ is given at high energy by 
the Casadei & Bindi fit and adding a DM contribution equal to (/)^+ + 0^+ . Solar modulation is implemented using 
the force field approximation [i^ with a Fisk potential of 300 MV. Fig. [T] illustrates the fit results in the case of 0=*= 
direct production. For PAMELA, both 100 GeV and 1 TeV WIMPs can accommodate the excess. As far as ATIC is 
concerned, we reproduce the observed feature in the case of a 1 TeV WIMP if we assume a DM clump with luminosity 
2.98 X 10^ Mq pc~'^ lying at a distance of 1.52 kpc from the Earth. The ATIC excess was reported in [3| where it 
was interpreted as evidence for a 620 GeV Kaluza-Klein species. We confirm that result in the case of a positronic 
line. In that case, no satisfying adjustment can be found adding a nearby subhalo, WIMP annihilations take place 
only inside a smooth Galactic DM distribution and the required cross-section <^ 10~^^ cm^ s~^, i.e. two orders of 
magnitude above our canonical value. As wc shall see in the following, the fit to the Fermi data points towards an 
incredibly bright clump. Indeed, the feature seen by Fermi is not very peaked and quite spread and requires, because 
of propagation effect, both a very large mass for the DM particle and a quite far away clump. We therefore exclude 
the combination of a thermal relic density plus a bright clump as a solution to this puzzle and we do not mention this 
case when dealing with other messengers. All parameters found in the best-fit cases are displayed in Tab. HI 
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TABLE I; Best fit values of the {D; L) couple in units of (kpc; Mq pc ^) for various DM particle masses and annihilation 
channels. 
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TABLE II: Upper rows : maximal values of the branching ratio Fp of WIMP annihilation into bb (100 GeV) or W^W 
(1 TeV) pairs allowed by the PAMELA antiproton data [S^]. Lower rows : the clump 7-ray flux above 0.1 GeV is expressed in 
units of 3 X 10"° 7 cm"^ s~^ [H]. 

DM annihilations within these clumps not only produce charged leptons but also antiprotons and 7-rays. We first 
allow a small branching ratio Fp into antiprotons either through the bb channel when the WIMP is light or through the 
W~'"W~ channel for a 1 TeV species. We compute the total antiproton flux (pp = (j)'^'^ + (l)p + <j>p with the same cosmic 
ray and DM models as for positrons. Requiring that the resulting signal docs not exceed the PAMELA antiproton 
data [2^ by more than Icr, wc derive the upper limits on Fp featured in Tab. |TT1 Antiprotons arc not forbidden if 
they are produced together with leptons but their abundance in the annihilation debris is sufficiently suppressed. 
We also present in Tab. |TT] conservative estimates for the detectability of the 7-ray emission from the lepton-fitting 
clumps of Tab. [D The 7-ray flux at the Earth is expressed in units of Fermi 5a sensitivity over 1 year of data taking 
for high-latitude point-like sources [2^. This is probably roughly close to the best the LAT instrument can do for 
a low-latitude source over its lifetime. We conclude that in all cases the LAT will detect the clumps. Notice that 
one of the leptophilic 1 TeV PAMELA cases has a very large flux of 1.11 x 10~^ 7 cm~^ s~^ and is already of the 
order of the EGRET point-source sensitivity of 2 x 10^^ 7 cm^^ s~^ [2^, and should have been already detected and 
correspond to one of EGRET unidentified sources. 

Fig.lHshows the probability of having the nearest DM clump of luminosity L within a distance D from the Sun. The 
abundance of nearby clumps and their properties are taken directly from the Via Lactea-II (VL-II) simulation [2^ . 
The high mass and force resolution, combined with a physical time step criterion [s^, allow VL-II to resolve subhalos 
even in the dense environment near the solar circle. The mean separation of subhalos with peak circular velocities 
Vmax > 5 kms""'^ is 9.6 kpc and their luminosities are 

L.mxion,|pc-»(^)'y^. (3) 
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FIG. 2: Best fit results in a clump luminosity-distance plane for different configurations, together with probabilities inferred 
from Via Lactea-II results. 

For comparison, the smooth VL-II main halo has a luminosity of L = 3.4 x 10^ Mq pc^"^, while the total luminosity 
is about 10 times higher (20| . At a given Vinax we assume a log-normal distribution of luminosities with factor of 3 
scatter, motivated by the substantial variance in the concentration cv found in nearby subhalos [20|. The bold line 
gives the median distances calculated from a random sample of observer positions. The long-dashed, dashed and 
doted lines stand for the 10th 1st and 0.1st percentiles, respectively. The points represent the locations of the best fits 
to the data in the L — D plane while the surrounding contours display the la excursions around these best fit values 
(as well as 3cr for ATIC and Fermi). We find that clumps fitting the PAMELA, ATIC and Fermi data arc far from 
the natural values indicated by VL-II. The most probable configuration is the PAMELA fit with a 100 GcV WIMP, 
which is inside the Via Lactca 3a contours. That configuration is found in 0.37% of all realizations. However, this 
scenario cannot accomodate ATIC data because is too small. Increasing the mass of the DM particle requires 
even brighter and less likely clumps at D ~ 1 kpc. As illustrated by the different PAMELA fit contours of Fig. [21 
the parameter degeneracy also increases as m-f^ gets higher. Basically, the PAMELA measurements do not constrain 
the spectral shape of the signal above 100 GeV and leave more lever-arm to the fits when WIMPs are heavy. For 
TeV WIMPs, there are clump properties which reproduce both the ATIC and PAMELA excess (see Fig 12]), and such 
a source would be well within the reach of Fermi (Tab. HII- However, in the standard CDM halo it is very unlikely 
to exist (p ~ 3 X 10~^). A DM spike or a higher cross-section would be required to get the needed luminosity from 
a smaller, more probable nearby subhalo.As shown in Tab. Ull the corresponding subhalo is within reach of Fermi. 
Concerning the best fit to the Fermi data, it points at the need for a clump that should be brighter than the whole 
Milky Way. Wc can safely associate a zero probability to this configuration. Note finally that the Via Lactca II 
contours extrapolate at lower values of the distance and clump luminosity. The corresponding clumps are not of 
particular interest for this particular study as for mean Via Lactea clumps, the natural luminosity decreases faster 
than what we gain from placing the clump closer. 

Our investigation of the Icpton anomalies in the presence of a nearby DM clump has finally led us to discover a 
subbtle interplay between the injection spectrum, the position of the substructure and the signal at the Earth. In 
particular, the commonly used criterion of a sharp cut-off as a smoking-gun signature of DM models producing Icptons 
(e.g. in Universal Extra Dimension Kaluza-Klein models [sH) is misleading if most of the DM signal stems from a 
highly clumpy halo. At TcV energies, leptons detected at the Earth arc produced within a distance of ~ 1 kpc and a 
clump becomes less visible if it lies outside that region. The opposite case is also possible since one can obtain very 
peaked Icpton spectra even with a soft injection spectrum, as long as the nearest clump is close enough. In general, 
as far as spectral shapes are concerned, it is important to bear in mind the spectral distortions induced by cosmic ray 
propagation. 

As an illustration, wc show in Fig. 121 how the signal from DM can exhibit a double-bump feature if, in addition to 
the contribution from a smooth DM halo, two nearby clumps are taken into account (solid light-colored line). In this 
example, the subhalos lie at a distance of 0.9 and 4.3 kpc from the Earth and their luminosities are of order 10* and 
10^" Mq pc^'^ respectively. Indeed, a specificity of lepton propagation is that regardless of the overall normalization. 
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FIG. 3: Electron-positron spectra resulting from one (dotted) or two (solid) nearby subhalos. 
DM species with thermal annihilation cross-section is considered. 



In the latter 620 GeV 



a feature observed at energy E can always be produced by a source at distance D with an injection energy Es as 
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However, without additional enhancements, having such bright nearby clumps 



is practically ruled out and we use them here for pedagogical purpose. CDM subhalos with L > 10**MqPC~'^ have 
VJnax > 25 kms^^ and host relatively bright dwarf galaxies [32| . If such dwarfs existed nearby, they should have 
been observed. As regards the light-colored dased curve of the right panel of Fig. [3l the sharp edge at 800 GeV is 
associated to a strong local DM annihilation and is produced by the VL-II smooth halo whereas the bump at ~ 100 
GeV comes from a single nearby clump located at 3.2 kpc. The cross section has been increased up to a value of 
10~^^ cm'^ s~^. This case seems somewhat more probable than the 2-clumps configuration. These examples illustrate 
how tricky boost factors are. Shifting upwards the DM cosmic ray fluxes turns out to be wrong especially in the light 
of the subbtle effects introduced by propagation. 

This letter is focused on the astrophysical boost associated to a nearby DM clump. We have restrained ourselves 
from using annihilation cross-sections in excess of the canonical value of 3 x 10^'^^ cm'^ s~^. As a final remark, we 
point out that all the fitted spectra would have been obtained with significantly less luminous subhalos should the 
annihilation cross-section be directly enhanced. If the Sommerfeld effect [33j_34| is at play, the signal from small 
clumps is enhanced with respect to the contribution from larger substructures 35[ since the velocity dispersion of DM 
particles decreases with the mass of the host subhalo. The blue contours of Fig. [2] are shifted towards smaller values 
of L and get nearer to the mean predictions of the VL-II simulation, with a much larger probability of occurrence (by 
simply assuming L ^ L x c/Vmax, the probability of ATIG best fit case increases from p ~ 3 x 10"^ to 14 percent!). 
Due to the lack of recent data, we had to do some assumptions on the electron spectrum and propagation parameters. 
Therefore new data may rule out some of these assumptions and hence lead to different conclusions. Future detailed 
measurements of the lepton spectrum from Fermi, and (on a much longer time-scale) of the positron fraction from 
e.g. PEBS [36| will provide us with much greater information on the nature and origin of the lepton anomalies, and 
perhaps shed light on galactic dark matter. 
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